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In magnetic superconductors vortices polarize spins nonuniformly and repolarize them when moving. At a 
low spin relaxation rate and at low bias currents vortices carrying magnetic polarization clouds become polaron- 
like and their velocities are determined by the effective drag coefficient which is significantly bigger than the 
Bardeen-Stephen (BS) one. As current increases, vortices release polarization clouds and the velocity as well 
as the voltage in the I-V characteristics jump to values corresponding to the BS drag coefficient at a critical 
current J c . The nonuniform components of the magnetic field and magnetization drop as velocity increases 
resulting in weaker polarization and discontinuous dynamic dissociation depinning transition. Experimentally 
the jump shows up as a depinning transition and the corresponding current at the jump is the depinning current. 
As current decreases, on the way back, vortices are retrapped by polarization clouds at the current J r < J c . As 
a result, polaronic effect suppresses dissipation and enhances critical current. Borocarbides (RE)Ni9B2C with a 
short penetration length and highly polarizable rare earth spins seem to be optimal systems for a detailed study 
of vortex polaron formation by measuring I-V characteristics. We propose also to use superconductor-magnet 
multilayer structure to study polaronic mechanism of pinning with the goal to achieve high critical currents. 
The magnetic layers should have large magnetic susceptibility to enhance the coupling between vortices and 
magnetization in magnetic layers while the relaxation of the magnetization should be slow. For Nb and proper 
magnet multilayer structure, we estimate the critical current density J c ~ 10 9 A/m 2 at magnetic field B » 1 T. 

PACS numbers: 74.25.Wx, 74.70.Dd, 74.25.Ha 



I. INTRODUCTION 



The conception of vortex as a polaron [ 1 1 was initiated by 
experimental data on the critical current in Er-borocarbide and 
first we discuss these data. The family of quaternary nickel 
borocarbides (RE)Ni2B2C (RE is rare earth magnetic ion) is 
an interesting class of crystals which exhibits both singlet su- 
perconductivity and magnetic order at low temperatures. (2J- 
4 1 A number of the crystals in that family develop antiferro- 
magnetic order below the Neel temperature 7#, which is be- 
low the superconducting critical temperature T c , Because the 
spatial periodicity of magnetic moments is well below the su- 
perconducting correlation length, superconductivity coexists 
quite peacefully with the antiferromagnetic order. In contrast, 
the ferromagnetic order, antagonistic to the Cooper pairing, 
leads to dramatic changes in both magnetic and superconduct- 
ing orders in the coexistence phase of singlet superconductors, 
for a review see Ref. Q. The compound ErNi2B2C with 
T c = 1 1 K and 7V = 6 K attracts numerous attention when 
it was realized that below the phase transition from incom- 
mensurate spin density wave (SDW) to commensurate SDW 
at T* = 2.3 K the phase with a weak ferromagnetic order- 
ing may emerge. |7) It was concluded that in ErNi2B2C 
below TV the incommensurate SDW develops with effective 
Ising spins oriented along the a-axis and with the wave vector 
Q = 0.5526 b* from neutron scattering measurements. (8l|9) 
Here b* = 2n/b and b is the lattice period along the b-axis. At 
T* the transition to the commensurate phase with Q = 0.55 b* 
leaves one out of 20 spins free of SDW molecular field. These 
Er spins with the magnetic moment yU = 7.8^g are easily po- 
larizable by the magnetic field along the a direction. The spin 
magnetization in the magnetic field H = 2000 G in tempera- 
ture interval 2 K - 4 K follows M sp /H « /jMJ{k B T), where 
M s * 56 G, see Fig. 4 in Ref. Q. This value, M s = fin, 



corresponds to magnetization when all "free" spins with the 
concentration n order ferromagnetically, The same value M s 
was obtained by extrapolation of the magnetization at temper- 
ature 2 K in fields H > 1500 G to H -> 0. d Nevertheless, 
the Hall probe measurements without an applied field below 
T* found internal magnetic field much lower than M s and no 
spontaneous vortex lattice was seen. IfTTI High polarizability 
of spin system in ErNi2B2C is a key point for our following 
discussion. 

As hope to observe remarkable consequences of weak fer- 
romagnetic phase coexisting with superconductivity waned, a 
puzzle on ErNi2B2C critical current behavior at low tempera- 
tures remained. It was discovered by measuring the hysteresis 
in the M - H loops and transport measurements that new pin- 
ning mechanism develops below 3 K for which the critical 
current increases as temperature lowers down to a: 1 .5 K fol- 
lowing approximately the enhancement of magnetic suscepti- 
bility. Qsma 

To explain these data the conception of vortex as a polaron 
was proposed, i.e. formation of polaron-like vortices dressed 
by the polarization cloud of magnetic moments. [ 1 1 Gener- 
ally, the polaronic mechanism is inherent to all magnetic su- 
perconductors but it is most pronounced when the magnetic 
system is highly polarizable, as in the case of ErNi2B2C below 
2.3 K. To clarify this mechanism, we recall that the magnetic 
field is nonuniform within the vortex lattice being strongest 
near the vortex cores. Consequently, the polarization of the 
magnetic moments is also nonuniform. When vortices move 
they should repolarize the magnetic system, otherwise they 
would lose the energy gained by polarization (the Zeeman en- 
ergy). The process of repolarization depends on the dynamics 
of magnetic system. In the following we consider the relax- 
ation dynamics of free spins in ErNi2B2C. The repolarization 
process is controlled by the relaxation time t which should be 
compared with the characteristic time a/v needed to shift the 
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FIG. 1 . (color online) Schematic view of the vortex lattice in the pres- 
ence of free Ising magnetic moments along the a axis. The vortex lat- 
tice is tilted from the applied magnetic fields in the ac plane due to the 
polarization of the magnetic moments. The vertical columns show 
the vortex cores. The polarized magnetic moments are nonuniform 
in space due to the spatial modulation of the vortex lattice magnetic 
field. Due to the Lorentz force F L vortices move along the x axis. In 
moving lattice, there is a phase shift between the magnetic induction 
(dashed line) associated with the vortex lattice and the magnetization 
(solid line) caused by the retardation in the response of magnetic mo- 
ments to the vortex magnetic field. 



vortex lattice moving with the velocity v by the vortex lattice 
period a - (<f>o/B) l/2 . Here (Do is the flux quantum, B is the 
magnetic induction and we assume a square vortex lattice. For 
t » a/v the magnetic moments slow down strongly the vor- 
tex motion. At some critical velocity and critical current J c , 
the vortices are stripped off polarization clouds. The corre- 
sponding jump in velocity is more evident for a large r's. As 
current decreases, the vortices become retrapped again at the 
current J r < J c . Since the voltage V cc v, the I-V character- 
istics show hysteresis. The physics here is similar to that of 
a polaron with vortices playing the role of electrons and the 
magnetic polarization in place of phonons IT3l . 



n. GENERAL EQUATIONS 

The ErNi2B2C crystals have orthorhombic structure below 
TV with domains where a- and fo-axes change by 90° in neigh- 
bouring domains. We consider a clean single-domain crystal. 
In multi-domain crystals, the domain walls also provide a pin- 
ning of vortices, which is sharply peaked when the vortex lines 
are aligned with the domain walls. [ 14 J We consider the vor- 
tex lattice induced by an applied magnetic field H tilted by the 
angle 6 with respect to the crystal c axis. As revealed by neu- 
tron scattering, vortices form square lattice in ErNi2B2C. iflTI 
We choose the z axis along the direction of vortex lines at rest 
and x axis in the ac-plane, see Fig. [T] The vortex line deviates 
from the applied field H due to the magnetic moments fl5l . 
The system is assumed to be uniform along the direction of 
vortex lines. In static situation the direction of vortex lines is 
determined by the effective field H + 47rM. Here M is the spa- 
tial average of the magnetization. We denote by a the angle 



between vortex lines and the c axis. 

The Lagrangian £{Rj(Y), M : (r, t)} for the whole system 
reads 



£{Ri(f), M(r, ?)} = £ M {M z (r, ?)} + £»{Hi(01 
+Xi„,{M z (r, t), R,(f)} + £„,{R,} + £ F {3}, 



(1) 



where X.m{M z (t, f)} is the Lagrangian for the magnetic sub- 
system 



£ M \M z (r, t)} = - j dr 2 M 2 (r, 



f)l(2 Xzz ), (2) 



and X v {Ri(0> r ji is the Lagrangian for the interaction between 
vortices and pinning potential due to quenched disorder 



(3) 



Here U (R, - r y) is the pinning potential at r ; . Further, _£ n ,{R,-} 
is the vortex- vortex interaction and Xf{J) = Z; J ' R^o/c is 
the Lagrangian due to Lorentz force in the presence of bias 
current density J. Here Xzz is the magnetic susceptibility at 
the working external magnetic field. It describes response of 
magnetic moments to nonuniform component of the field in- 
duced by vortices. 

In the London approximation the magnetic field of the vor- 
tex lattice inside the crystal is (r = x, y) 



B z (r) = B,^ 



cos(G • r) 
A 2 G 2 + 1 ' 



(4) 



where B z is the averaged magnetic induction, G are recipro- 
cal vectors of the square lattice and A is the superconducting 
penetration length renormalized by the magnetic moments. It 
is given by the expression A 2 = A 2 L (l - Anx zz ), where A L de- 
scribes magnetic field penetration in the absence of the mag- 
netic moments. Il5l [T6l4T9ll Note that the magnetic susceptibil- 
ity^ = M z /B z is smaller than 1/(4tt), i.e. Xa < 1/(4tt). The 
magnetic fluctuations (M,M Z ) ~ Xzz/i^ ~ At7l Xzd diverge when 
Xzz — * 1/(47!"), which indicates instability of the magnetic sys- 
tem ll20ll . Here we also ignore anisotropy of the penetration 
length. 

In the Lagrangian the interaction between vortex line at 
R, = (Xj,yi), and the magnetic moments is determined by the 
term 



Xi„,{R,-,M 



dr 2 B z (Ri - r, t)M z (r, t), (5) 



where we describe the magnetic moments in the continuous 
approximation via the magnetization M z (r, t), because dis- 
tance between free spins, 35 A, [9 1 is much smaller than the 
London penetration length A, about 500 A. ifTBI We ignore the 
pair breaking effect of the magnetic moments because they 
suppress Cooper pairing uniformly as distance between free 
spins is much smaller than the coherence length, and thus the 
pair breaking effect by the moments does not introduce pin- 
ning. 
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For 



„ - the system shows hysteresis in the I-V curve while for 
F p = 2, no hysteresis is present. The dark yellow dotted line denotes 
the unstable solution. 



Both the magnetization and vortices are governed by a re- 
laxation dynamics characterized by the dissipation function 
7?{R/(f), M z (r, t)) = <R M . + n v , where 



RM-\M z (r)} 



hi 



dr 2 M 2 (r), 



(6) 



Here t is the relaxation time for a single spin and 77 is the 
Bardeen-Stephen drag coefficient per unit vortex length, 77 = 
Oy/(27r^ 2 c 2 p„) with p n the normal resistivity slightly above 
T c . The equation of motion for vortices is given by the Euler- 
Lagrangian equation of motion 



d S£ _ 8Z + 8R _ 
dtSRi SRi 5Ri ~ 



(7) 



which gives 



<9R; 5X VV {R ; ,R ; } 3i:i nt {R ; ,M} 
] ~df = " 



dRi dRi 
dU(R t - rj) 



■S 



dR t 



+ F 7 



(8) 



with Fl = OoJ/c being the Lorentz force. 

We neglect here the effect of quenched disorder because the 
vortex motion quickly averages out the disorder and the lattice 
ordering is improved ET1I221 . In the lattice phase _£, T = due 
to symmetry. The equation of motion for vortex lines then is 



dRt d£URi,M z 



+ F, 



' dt dRi 
The dynamics of the magnetization is governed by 



dM : (r, t) 
3t 



M,(r, t) 



B z (r) 



(9) 



(10) 



From Eq. ( 10 1 we see that relaxation time of the magnetization 
measured experimentally in the crystal is Xzz T - The force due 
to magnetic moments is the same for all lines and the vortex 
lattice moves as a whole. The motion of vortex lattice center 
of mass, w(f), along the x-axis is described by the equation 



du 
7 di 



d_ 
du 



drB z (x + u,y, t)M z {r, t) 



+ F, 



(11) 



Using the linear response approach to relate magnetization 
with the magnetic field, we obtain 



'1 



du d 
dt du 
B z (r',t') + F L 



J*drdr'B z (x + u,y, t) J dt'x. 



zz (r -r',t- f) x 
(12) 



The vortex lattice moves with a constant velocity, u = vt, in 
the steady state t » r. Integrating over coordinates and time 
we obtain 



rjv • 



z 



(A 2 G 2 + l) 2 



+ Ft 



(13) 



where ^ zz (k, u) is the dynamic magnetic susceptibility in the 
Fourier representation. We see that the magnetic moments 
affect strongly vortex motion if a) the resonance Cherenkov 
radiation condition v ■ G = Q(G) is fulfilled, where O(k) is 
the frequency of magnetic excitations with the momentum k 
and £2(k) » F(k), where F(k) is the relaxation rate of ex- 
citation, and b) dynamics of magnetic system is dominated 
by relaxation, O(k) < F(k). In the former case, discussed in 
Refs. [23—25], the magnetic moments renormalize the vortex 
viscosity at high velocities when alternating magnetic field of 
vortices is able to excite magnons. Here we consider the latter 
case of free moments described by the relaxation dynamics 
according to Eq. (TTTJb with^ z ,(k, a>) given by 



Xzz(k,o)) =^-sin a 



1 



X = 



k R T 



(14) 
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FIG. 3. Dependence of the critical current J c and retrapping current 
J r , and corresponding electric fields E c and E r on F p . When F p < 8, 
hysteresis in I-V curve disappears. 
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at temperatures T below 3 K for ErNi2B2C. 

Renormalizing time in units of t\, length in unit of 1/Gi, 
force per unit vortex length in unit of t]I{tG\x) we obtain 
equation for velocity 



v + F 



P„2 



V 2 + 1 

where we have accounted only for dominant lattice wave vec- 
tor Gi = (2n/a,0,0) and introduced magnetic pinning force 
per unit vortex length F p = ti^TX 2 sin 2 a/(4-7i 2 A 4 ri). At low 
bias current (low F£) velocity is proportional to Fi but with 
enhanced effective viscosity, v * Fl/(1+F p ). At a large Fi (a 
large v) renormalization vanishes and I-V characteristics be- 
comes the usual Bardeen-Stephen one. Importantly, at F p > 8 
change occurs through sharp transition, as shown in Fig. [2] 



Equation ( 15 1 at F„ > 8 for an intermediate J has three real 



solutions, the highest V3 corresponds to decoupled motion of 
vortex lattice and magnetization, the lowest vi corresponds to 
the motion of vortex-polaron, and the intermediate solution V2 
corresponds to an unstable state. 

The jump at J c , identifying experimentally as depinning 
transition, is caused by the dissociation of the vortex-magnon 
polaron. It is very similar to the dissociation of usual electron- 
phonon polaron in high electric fields as described theoreti- 
cally |26| and confirmed experimentally in metal oxides IT271 . 
Upon decreasing the current the vortices are retrapped by the 
polarization clouds at a threshold current J r and the vortex 
lattice moves with a significantly enhanced viscosity at lower 
currents. The calculated J c , J r and corresponding electric 
fields are shown in Fig. [3] At large F p the critical current 
is independent of r. 



0.03 



Xc<S>o sin a 
G1A 4 



J c decreases with temperature as J c ~ l/T and decreases with 
the magnetic field as J c ~ 1 / VB. 

Let us explain the origin of the jumps at J c and J r . The 
dependence of the magnetization on the velocity of moving 
vortices is 



length £ « 13 nm [ 15 1 and the normal resistivity p„ = 5 fiQ.-cm 
at T c . [28 1 The relaxation time x T m ErNi2B2C is long be- 
cause the dynamics of majority of spins is strongly suppressed 
by the formation of the SDW molecular field as was found by 
the Mossbauer measurements.! 29 1 The relaxation time drops 
(15) very fast below 10 K and reaches the value \ T ~ 5 x 10~ 10 s 
at T — 5 K. The data at lower temperatures however were 
not measured. Thus the only information we have so far is 
F p > 50 sin 2 a. 

The critical current for ErNi2B2C reported in Ref. IfTUl is 
about 250 A/cm 2 for B — 0.1 T, T — 2 K. which corresponds 
to a — 2.5° according to Eq. ( [To} . The applied magnetic 
field was close to the c axis in experiment, but the precise 
angle 8 was not reported. flTOll . The estimate of the order 1° 
is reasonable, but the quantitative comparison is not convinc- 
ing as we do not know r and thus F p below 2.3 K. We pre- 
dict hysteretic behavior in ErNi2B2C, strong dependence of 
voltage and of the critical current on the angle 8, at least for 
8 » 0.15°. Hence, the real check of polaronic mechanism 
should be by measuring the I-V characteristics. We estimate 
the critical current reaches values as high as 10 6 A/cm 2 at high 
angles at T = 1 K and B = 0.1 T. 

The effect of ordered spins on the vortex motion is similar 
to that described in Refs. 23 -25 for an antiferromagnet. When 
the Cherenkov radiation condition v ■ G = Q(G) is satisfied, 
excitation of magnons results in enhanced drag coefficient by 
transferring energy from vortex motion to the magnetic sub- 
system. This occurs at high velocities (high currents) of vor- 
tices, due to a gap in the magnon spectrum and a large velocity 
of magnon, leading to a voltage drop in comparison with the 
bare BS result. 

In the incommensurate SDW when T > T* , some spins 
experience quite weak SDW molecular field. Thus, they are 
polarized by vortices and exhibit polaronic effect and pinning. 
This accounts for the increase of pinning in ErNi2B2C as T 
decreases below T N , see Ref. IfTOll , and also the pinning in the 
holmium borocarbide below 7V QUI 



(16) 



M,(r, v, t) = xB sin 2 a } i 



cos[G-r-/3(v)] 



01 2 G 2 + 1)[1 +{G x vtx) 2 ? 11 



(17) 

with tan(y6) = G\vrx- Nonuniform component of the magne- 
tization and thus the polarization effect decrease with veloc- 
ity. On the other hand, the retardation between the magnetic 
field and the magnetization, as described by the phase shift 
B(v), increases with the velocity. This positive feedback and 
the increase of retardation with velocity ensure discontinuous 
transitions at J c and J r . 



HI. DISCUSSION OF EXPERIMENTAL DATA FOR 
ERBIUM BOROCARBIDE 



A large parameter F p is required to have a strong pinning 
due to the polaron mechanism. It is expressed in terms of t 
as F p w 10 n ^rsin 2 as -1 , where we have used the coherence 



Next we discuss the effect of quenched disorder. In the 
presence of quenched disorder, the vortex lines adjust them- 
selves to take the advantage of the pinning potential, which 
destroys the long-range lattice order. Below a threshold cur- 
rent, vortices remain pinned (actually they creep between pin- 
ning centers due to fluctuations), where the polaronic mech- 
anism does not play a role in this region. When the current 
is high enough to depin the vortices from quenched disorder, 
vortices start to move and the lattice ordering is improved. 
The vortex viscosity is enhanced by formation of polaron with 
the nonuniformly induced magnetization. The polaron disso- 
ciates and the system jumps to the conventional BS branch 
at a critical velocity (current). Pinning due to quenched dis- 
order works in the static region and polaronic pinning works 
in the dynamic region. The critical current of the whole sys- 
tem therefore is the sum of these two threshold currents. Note 
that magnetostriction in combination with quenched disorder 
enhance the polaronic pinning mechanism. 
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IV. RESPONSE OF THE VORTEX LATTICE TO AN AC 
DRIVING CURRENT 

Here we study the response of a vortex-polaron to an ac 
driving current [31 1. We write the equations of motion for 
magnetization m(f) - M(G\,f)A 2 G 2 l$>QX an d me vortex lat- 
tice center of mass u(t)) as 

d,m(t) = ~[m(t) - exp[-;'w(f)]], (18) 



d,u = F L - Im [F p exp(/«)m(f)] , (19) 



with an ac Lorentz force Fi 
we obtain equation for u(t): 



F dc sm(a>t). Eliminating m(t) 



du 
dt 



> f dt' 
Jo 



sin[u(t) - u(t')] exp(f' - t). 



(20) 



First we consider an ac current regime with a low amplitude 
^ae/[w(l + F p )] <K 1. Then the vortex lattice oscillates, u = 
Re[u ac exp(iajt)], with the amplitude 



U ac = F ac0'?7effW + a p) 

?7eff = l+F p (u 2 + l)- ] 



(21) 

a p = F p aj 2 (aj 2 + (22) 



For a high frequency cj » 1, the effect of magnetization is 
to introduce the pinning potential Um = F p u 2 /2 with strength 
F p . In this case, the vortex lattice follows the driving force 
much faster than magnetization, which remains almost time 
independent. The polarization of the magnetization results in 
periodic pinning potential with the periodicity of vortex lattice 
as it was induced by the same lattice at previous positions and 
previous moments of time. For a low frequency a> «: 1, the 
effect of magnetization is to renormalize the drag coefficient 
from r] to ?7eff = 1 + F p . In this polaron region, the magneti- 
zation follows vortex motion by formation of vortex polaron, 
as in the dc case to — 0, resulting in enhancement of viscosity 
and suppression of ac dissipation. The dissipation power of 
the whole system, averaged over time, D(co) = (Fi(t)v(t)),, 
is reduced due to the presence of magnetic moments. In the 
linear region with vortex polaron, we obtain 



D(co) = 



W 2 7? e ff 



(23) 



This dissipation power should be compared with the case 
without magnetic moments (at F p = 0), Do = F\j2. For 
w » 1 we get D/Dq * 1 and for co <k 1 we get D/D {) 
( I ■ Fp)~ l . The frequency dependence of the normalized dis- 
sipation power D(oj)/D(), effective viscosity 77^, and pinning 
strength a p is shown in Fig. |4] The dissipation of the system 
in the presence of magnetic subsystem is strongly reduced in 
linear regime F ac < Fi c , which might be useful for applica- 
tions. 

Next we consider stronger driven force amplitude. In this 
hysteretic regime, we describe the system analytically in the 
adiabatic limit, w « 1. At the time moment t c , when F L (t) - 
Flc ~ 0.5F p , polaron dissociation leaves magnetization and 
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FIG. 4. (color online) (a) Dependence of the normalized dissipation 
power, D((±))/D , and (b) the effective viscosity rj cSf , pinning strength 
a p on the driving frequency u> in the linear regime F ac <k a>( \ + F p ). 
Here F„ = 20. 



vortex lattice weakly coupled because lattice moves now with 
a high velocity. The magnetization component m(t) after that 
moment relaxes as m( f) - exp(-f + f c ), and motion of vortex 
lattice is determined by the equation 



du 
~di 



= Fic + F p sin(t<) exp(-f + t c ). 



(24) 



When t — t c < 1 the velocity of the vortex lattice oscillates 
with the frequency £2 = F ic, 



v = F Lc + F p sm(Fu;t) exp(-f + t c ), 



(25) 



but oscillations relax on the time scale of unity. These post 
dissociation oscillations are caused by the motion of vortex 




-20 -15 -10 



-5 5 

F L [Tl/CGjtx)] 



FIG. 5. (color online) Dependence of the vortex velocity v(f) on the 
driving force F L (t) = F ac sin(wr) with F ac = 20. Here F p = 20. 
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FIG. 6. (color online) Time evolution of the vortex velocity, v(f), and 
of magnetization \m{t)\ in the presence of ac driving force F L {t) = 
F ac sin(wf) at several frequencies to. (a): to = 10, (b): to = 1, (c): 
w = 0.1. We take F ac = 20 and F„ = 20. 



lattice in periodic potential induced by remnant retarded mag- 
netization. 

To take both the retardation and nonlinearity into account 
for an arbitrary to we solve numerically Eqs. ( fl~8] > and ( 19 1. 

> 



We consider the interesting region F p > 8, where the disso- 
ciation of vortex polaron is possible due to nonlinear effects 
at u > I. We take F p = 20 in the following discussion. The 
hysteretic behavior of vortex lattice velocity vs. driving force 
is shown in Fig. [3] At frequencies to < 1 which are simi- 
lar to the dc case to = 0, we see the following sequence of 
events during the period of Fi(t): polaron formation near low 
\F L \ (interval of low vortex velocity), then polaron dissociation 
(velocity sharp increase) followed by region of vortex oscilla- 
tions on the background of average high velocity, decrease 
of velocity as the Lorentz force drops and vortex retrapping 
(sharp drop in vortex velocity) and then dissociation again at 
negative -Fi c (sharp drop in velocity). The results for behav- 
ior of vortex velocity in time, v{t), at F ac — 20 > Fi c and 
different to are shown in Fig. [6] 

At all frequencies to < 1 we see post dissociation oscil- 
lations caused by motion of decoupled vortices with respect 
to periodic potential created by nonuniform magnetization in- 
duced by the same lattice just before decoupling (when ve- 
locity was still low) and frozen for some periods of time af- 
ter decoupling due to the retardation effect. This self-induced 
pinning due to the retardation, and the amplitude of corre- 
sponding vortex oscillations reach maximum at a> « 1 . In a 



rough approximation we describe them by the equation 

du 
dt 



Flc + F p m d sin(w - u d ), 



(26) 



assuming approximately constant m and Fi = Fi c in the re- 
gions of maxima and minima of the Lorentz force. Here u d 
and m c { are the position of vortex lattice and the amplitude of 
magnetization at the moment of decoupling. This gives ap- 
proximate solution: 



v(t) ~ F Lc + F p m d sin(F Lc f), 



(27) 



which provides rough estimate for the oscillation frequency, 
£2 ~ Fic, when number of oscillations of velocity per the half 
period of Fi£t) is significantly bigger than unity. This expres- 
sion for the frequency in original units reads Q » InF^/aq. 
Such relation is anticipated for decoupled vortex moving in 
the pinning potential with periodicity a. 



V. ENHANCEMENT OF CRITICAL CURRENT DENSITY 
IN SUPERCONDUCTING/MAGNETIC MULTI-LAYERS 

The polaronic mechanism of pinning is promising for 
achieving a high critical current. We propose to use a su- 
perconducting (S) and magnetic (M) multilayer structure as 
shown in Fig. |7]to optimize such pinning mechanism. B2l To 
achieve high critical current the magnetic layers should have a 
slow relaxation of the magnetization. Secondly, the magnetic 
layers should have a high magnetic susceptibility at working 
magnetic field to ensure a strong coupling between magnetic 
moments and vortices. Thirdly, the penetration depth of the 
superconducting layers should be small, such that the magne- 
tization polarization varies rapidly in space. 

The vortex lattice is induced inside the S layers under ex- 
ternal magnetic fields. The vortex lattice moves in response 
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FIG. 7. (color online) Schematic view of multi-layer structure con- 
sisting of alternating magnetic (M) layers (purple) with thickness d„, 
and superconducting (S) layers (blue) with thickness d s . 
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to the Lorentz force when a transport current is present. In 
the quasistatic approximation, the motion of vortex lattice is 
given by 



l 2 VxVxB+B = (Do^^r- r,(f)] z, 



(28) 



where z is the unit vector along the z axis and r,(f) = ro - vf is 
the vortex i coordinate. The magnetic field inside the M layers 
is governed by the Maxwell equations 



V x (B - 4ttM) = 0, V • B = 0. 



(29) 



The dependence of magnetization M on B is determined by 
the material properties. With a strong field and in static 
case, M is a nonlinear function of B and generally can be 
expressed as M(r) = Jc/r 3 /[r - r',B(r')]- The character- 
istic length of magnetic subsystem is much smaller than A 
and we use a local approximation /[r - r',B(r')] = <5(r - 
r')/(B(r')- Induction B(r) has component uniform in space, 
Bo, and the nonuniform component nonuniform, B(r) <k Bq. 
Thus the spatially nonuniform magnetization M(r) is M(r) at 
[d/(B )/dfio]B(r) = # (B )B(r). In the following we con- 
sider isotropic magnetic subsystem characterized by a suscep- 
tibility ^o(Bo) at Bo in static case. The magnetic field inside 
the M layer is determined by the equation V 2 B = 0. Since 
only the spatially nonuniform component M and B are re- 
sponsible for pinning, we will focus on the nonuniform com- 
ponents in the following calculations and omit mark tilde. At 
the interface between the M and S layers, we use the standard 
boundary condition for the field parallel to the z-axis B z and 
field parallel to the interface Z?" 

B : \ s = B Z \ M , B"| s = (1 - 4 W) )B"|m. (30) 
Then we obtain the magnetic field inside the M layers 

mtn^n \ r -Gfe-+<4)1 exp (-iG x v x t) 
B m (G >Q,z) = a [e uz + e ° u *" ; J - ^ , (31) 



BUG>0,z) 



- e 



-G(z'+d,„) 



O exp(-iG x v x t) 



1 + A 2 G 2 



(32) 



e d m G(_ l+e d s kj x , 



(1 - X ')(e d ° k ° - e Gd -) + (1 + x')(l ~ e d "> G+d - k -y 



with z' - z - n(d s + d m ), k s = V/L 2 + G 2 , and^' = (1 - 
4nxo) ks/G. Here n is the layer index and the vortex motion 
is assumed to be along the x direction. We consider square 
lattice G = (m x 27t/a,my2n/a) with a = V^o7§o the lattice 
constant and m x , m y integers. 

We assume a relaxational dynamics for the M layers, 
M(a>) = ^(w)B m («), with a dynamic susceptibility governed 
by a single relaxation time xot as in Eq. (jT0]l. In the steady 
state, we have 



M(G 



Jo 



exp[(f' - f)/Cror)]B m (G,z, t')dt'. (33) 



M depends on the history of vortex motion for a slow relax- 
ation of magnetization, and there is retardation between the 
time variation of induced nonuniform magnetization and vor- 
tex motion. As a result, the magnetization exerts a drag force 
to the vortex which is opposite to the driving force. The pin- 
ning force acting on a single vortex due to the induced magne- 
tization in one M layer is given by Fm - d ro J dxdy J dzM ■ 
B,„, which yields 



2a 2 xo GvxqtQI 



(\+A 2 G 2 fa 2 1+(Gvx t) 2 ' 
(34) 

The TV curve is determined by the equation of motion for 
vortex d s T]v = d s Fi - Fm with the electric field E = Bv/c and 
the Lorentz force F L = J<f>o/c. We consider a realistic case 
where a/(2n) <K d m , d s . Taking into account only the domi- 
nant contribution G x = 2n/a and G v = in the summation, 
we obtain the same equation as Eq. (IB) , but with a different 
parameter 

f - — I 1 \ 2 *t!% 

" 2r,d s \l-2nxo) A\2ny ^' 



after introducing the same dimensionless units as before. 

Hysteresis is developed when F p > 8. For typical parame- 
ters for Nb superconductor £ * A as 40 nm, p n as 10 -6 Q. ■ m 
and a = 40 nm at B as IT and xo - 0.05, F p > 8 requires 
Xot > 1 ps. For the relaxation time of order xqt ~ 1 y"s, the 
effective viscosity is enhanced by a factor of 10 6 compared to 
the bare BS one at v < a/Ct'oT). The effective critical current 
density for the whole system is given by 



1 



(36) 



For d s = d„, 
retrapping current J r is 



2-4nx a j (2n) 4 A 4 d s + d n 
100 nm , we obtain J c * 10 9 A/m 2 . The 



Jr 



1 



rjad s 



1 



1 - 2nxo V tit A 2 4n 2 d s + d m 



(37) 



For the parameters used before and xo T = 1 M s > we estimate 
J r as 2 x 10 6 A/m 2 . 

Let us discuss the optimal materials for the S and M lay- 
ers. Superconductors with a smaller A are preferred because 
the critical current decreases as A~ A . The smaller A, the more 
nonuniform is the magnetic field distribution inside the M lay- 
ers, hence stronger pinning. The viscosity in the branch with 
vortex polaron is proportional to r while the critical current 
does not depend on t for sufficiently large r. The slow mag- 
netic dynamics can be realized in certain spin glasses, where 
the relaxation of magnetization is governed by a broad spec- 
trum of time scale, with average time of the order 0.1 fis 
Glllll. For CuMno.os, xo ~ 0.002 atB= 1 T. (35j One may 
enhance xo by tuning the concentration of magnetic metal in 
alloys. [36 1 One may also use superparamagnets with r as 
large as 1 s and with a huge xo due to large magnetic mo- 
ments in superparamagnets [ 37 -39 ] and the recently synthe- 
sized cobalt-based and rare-earth-based single chain magnets 
with xo ~ 0.05 at /J = IT and 10~ 6 s < xqt < 10~ 4 s. R0H4H 
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Now we discuss the optimal thickness of M and S layers. 
For d m » a, we have B m (G > 0) ~ exp(-2nd m /a) when 
—d m <K z' <sc according to Eqs. (31 1 and (32i. The mag- 



netic induction and the magnetization are almost uniform in 
the lateral direction in the middle of the M layer. As a result, 
the pinning force becomes practically d m independent in this 
case. In other words, the pinning is effective only near the 
boundaries between S and M layers in the area of thickness 
of the order a. On the other hand, the Lorentz force is pro- 
portional d s . Thus the effective critical current of the whole 
system J c is proportional to l/(d s + d m ) as described by Eq. 



(36 1. Therefore the thinner of both M and S layers, the higher 



is the critical current of the system. 

The M/S multilayer structure is naturally realized in cer- 
tain superconducting single crystals, such as (RE)Ba2Cu3C>7 
EH S3 and RuSr 2 GdCu 2 8 06). For (RE)Ba 2 Cu 3 7 , mag- 
netic RE ions interact weakly with superconducting electrons 
because they are positioned between the superconducting lay- 
ers. They order at very low Neel temperatures of the order 
TV ~ 1 K. The polaronic mechanism is important at T > T^, 
where spins are free. The London penetration depth of cuprate 
superconductors is large A m 200 nm, thus the critical cur- 
rent is reduced significantly compared to that for Nb multi- 
layer structure, because J c drops as I /A 4 . Another natural 
realization is the recently discovered iron-based superconduc- 
tors, such as (RE)FeAsOi_ x F x , where RE ions ordered anti- 
ferromagnetically below TV ~ IK. BTll In RuS^GdCuaOg 
the magnetic moments order ferromagnetic ally above T c thus 
the dominant enhancement of vortex viscosity is due to the 
radiation of magnons IT23TJ251 . 



VI. CONCLUSIONS 



polarized magnetic moments at the retrapping current which 
is smaller than the critical one. The results of such polaronic 
mechanism are in qualitative agreement with the experimen- 
tal data IfTUl IT31 but measurements of the I-V characteristics 
are needed to establish the quantitative agreement and confirm 
validity of such a model for Er borocarbide. The polaronic 
mechanism should be at play also in Gd and Tb borocarbides 
superconductors in the commensurate SDW phase and strong 
effect may be present in Tm borocarbide above 7V 

We derive the response of the magnetic superconductors in 
the vortex state to the ac Lorentz force, F L {t) — F ac sin(wf), 
taking into account the polaronic effect. At low amplitudes 
of the driving force F ac the dissipation in the system is sup- 
pressed due to the enhancement of the effective viscosity at 
low frequencies and due to formation of the magnetic pinning 
at high frequencies a>. In the adiabatic limit with low frequen- 
cies oj and high amplitude of the driving force F ac , the vortex 
and magnetic polarization form a vortex polaron when F i(t) 
is small. When Fi increases, the vortex polaron accelerates 
and at a threshold driving force it dissociates, i.e. the vortex 
motion and the magnetization relaxation decouple. When Fi 
decreases, the vortex is retrapped by the background of rem- 
nant magnetization and they again form vortex polaron. This 
process repeats when Fi{t) increases in the opposite direc- 
tion. Remarkably, after dissociation, decoupled vortices move 
in the periodic potential induced by magnetization which re- 
mains for some periods of time due to retardation of mag- 
netization after the decoupling. At this stage vortices oscil- 
late with high frequencies determined by the amplitude of the 
Lorentz force at the moment of dissociation. 

We propose to make multilayer system M/S where one can 
optimize superconducting and magnetic layers to achieve high 
critical current. 



In conclusions, vortices in magnetic superconductors polar- 
ize magnetic moments and become dressed and polaron-like. 
At low currents and long spin relaxation time the nonuniform 
polarization induced by vortices slows their motion at currents 
for which pinning by crystal lattice disorder becomes ineffec- 
tive. As current increases above the critical one, vortices re- 
lease nonuniform part of the polarization and the velocity as 
well as the voltage in the I-V characteristics jump to much 
higher values. At decreasing current vortices are retrapped by 
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